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ABSTRACT. We study the dyadic model of the Navier-Stokes equations introduced by 
Katz and Pavlovic. They showed a finite time blow-up in the case where the dissipation 
degree a is less than 1/4. In this paper we prove the existence of weak solutions for all a, 
energy inequality for every weak solution with nonnegative initial data starting from any 
time, local regularity for a > 1/3, and global regularity for a > 1/2. In addition, we 
prove a finite time blow-up in the case where a < 1/3. It is remarkable that the model 
with a = 1/3 enjoys the same estimates on the nonlinear term as the 4D Navier-Stokes 
equations. Finally, we discuss a weak global attractor, which coincides with a maximal 
bounded invariant set for all a and becomes a strong global attractor for a > 1/2. 



1. Introduction 

The regularity of the 3D incompressible Navier-Stokes equations (NSE) remains a sig- 
nificant problem. This, among many other open problems connected with the 3D NSE, 
depends on the estimates on the inertial term (u • V)u in the equations. In this paper we 
study a dyadic model, which has similar properties to the 3D NSE, the same estimates on 
the inertial term, and the same open question concerning the regularity of the solutions. 

There have been many simple models proposed in the literature that capture some es- 
sential features of the 3D NSE. Among these are shell models of turbulence, which have 
been investigated for many years (see l2lllllfT2lll6M18in . Recently, some of these models, 
as well as some new ones, were extensively studied analytically. In |7|, Constantin, Lev- 
ant, and Titi study the "sabra" shell model of turbulence, proving a global regularity and 
the existence of a finite dimensional global attractor and inertial manifold. 

In 1141 . Katz and Pavlovic introduced another shell-type model, the dyadic model for 
the Euler and Navier-Stokes equations. This model, motivated by 1131 . is an infinite system 
of nonlinear ODEs that describes evolutions of wavelet coefficients. In 1 9 1, Friedlander and 
Pavlovic proposed a three-dimensional vector model for the Euler equations, similar to a 
quasi-linear approximation of the 3D Navier-Stokes system constructed by Dinaburg and 
Sinai 1 8 1. Both of these dyadic models can be reduced to the following system of nonlinear 
ODEs: 

(1.1) ^-u n + v\ 2an u n -\ n ul^ + \ n+1 u n u n+1 =g n) neN, 

dt 

where uq = 0. Here, A > 1, v > is the viscosity, and a > is the dissipation degree. 
Note that we also include a force g = [g\, g2, . ■ ■ ) in the model. For u = (iti, «2, ■ ■ ■), the 
dyadic model can be written as 

— u + vAu + B(u, u) = g, 
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where 

(Au) n = \ 2an u n , (B(u, u)) n = -\ n u 2 n _ 1 + X n+1 u n u n+1 , 
and Mo = 0. Note that for a = 2/5 the following are sharp estimates on the inertial term: 

\{B(u,u),Au)\ < lAul 3 / 2 ^ 1 / 2 !*! 3 / 2 , 

where (■, ■) and | ■ | are the Z 2 -inner product and norm respectively. The best known esti- 
mates on the inertial term of the 3D NSE are the same, with (•, •) and | • | being the L 2 -inner 
product and norm. 

Katz and Pavlovic [ 14 1 proved that under certain assumptions on the initial conditions, 
solutions of the inviscid dyadic model blow up in finite time in a norm stronger than the 
/ 2 -norm. Later, Waleffe |21 1 derived the inviscid dyadic model from the Burgers equation 
reducing it to ( II. II with v = 0. Recently, Kiselev and Zlatos 1 15 1 sharpened the blow-up 
result for the dyadic model and studied a very similar Obuhov model (see 1 17 1) proving a 
global regularity of every solution with regular initial data. In addition, the existence of a 
global attractor of the inviscid dyadic model is proved in 1 5 1 . This surprising fact is a result 
of a self-dissipation mechanism due to the loss of regularity of solutions. 

In this paper we will study the viscous dyadic model, i.e., the model ( II. II with v > 0. 
In this case Katz and Pavlovic 1 14 1 obtained a finite time blow-up of solutions with certain 
initial data when a < 1/4. Our main goal is to prove a finite time blow-up in the case 
where a < 1/3. It is remarkable that in the critical case a = 1/3 the following are sharp 
estimates on the nonlinear term: 

\(B(u,u),Au)\ < \Au\ 2 \A 1/2 u\. 

Note that the best known estimates on the nonlinear term of the 4D NSE are exactly the 
same. 

This paper is structured as follows. We start with surveying the dyadic model in Sec- 
tion|2] In Section[5]we introduce a functional setting and define weak and strong solutions. 
In Section|4]we derive some a priori estimates and prove the existence of weak solutions 
to the dyadic model. This is done by taking a limit of the Galerkin approximation, which 
also results in the energy inequality for a limit solution (which might not be unique) for 
almost all time. Then we show that the lack of backward energy transfer implies that every 
weak solution with nonnegative initial data satisfies the energy inequality starting from any 
time. Finally, using a classical NSE technique, we show a local regularity for a > 1/3 and 
a global regularity for a > 1/2. 

In Section [5] inverting the Sobolev-type estimates, we prove that every solution with 
large i/ e -norm blows up in finite time in i/ 1 / 3+e -norm, e > 0. Here, i/ 7 -norm is defined 
as \\u || 7 = (EA 27 "w 2 ) 1/2 . Note that we also use such a technique in |4|, where we study 
a similar model that has coefficients growing as power functions. That model is introduced 
as an example of a NSE-like dynamical system that possesses a weak global attractor, on 
which all the solutions blow up in finite time. The reason for the power-law growth for 
the coefficients is to mimic the growth of eigenvalues of the Stokes operator in 3D. It is 
remarkable that due to a slower growth of the coefficients, the model in |4| still possesses 
a gap between the regions of a local regularity and finite time blow-up. 

Lastly, in Section[6]we discuss a weak global attractor for the dyadic model. The weak 
global attractor is the minimal weakly closed weakly attracting set. Using results from 
0EL we show that the weak global attractor is also the maximal bounded invariant set. 
Moreover, for a > 1/2 all the trajectories are continuous in I 2 , which implies that the weak 
global attractor is in fact a strongly compact strong global attractor. 
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Note that there is still a gap between the regions of global regularity and blow-up in 
finite time, which means that the developed technique is not sharp enough to separate these 
two behaviors. Since most of the proofs in the theory of the Navier-Stokes equations go 
through for the dyadic model, a better understanding of the dyadic or similar shell models 
might provide insight into the regularity problem for the Navier-Stokes equations. 

2. Dyadic model 

Here we recall a derivation of the dyadic model for the equations of fluid motion by Katz 
and Pavlovic 1 14|. A cube Q C K 3 is called dyadic if its side length is 2 l , and the corners 
are on the lattice 2 Z Z 3 , for some integer I, For a dyadic cube Q with side length 2~' J , its 
parent Q is a unique dyadic cube with side length 2~ J+1 that contains Q. For m > 1, let 
C m (Q) be the set of all mth order grandchildren of Q, i.e., all the dyadic cubes with side 
length 2~i~ m that are contained in Q. For instance, C 1 (Q) consists of 2 3 children of Q. 

Now a scalar-valued function u(x, t) can be represented by the following wavelet ex- 
pansion 

u(x,t) = ^2u Q (t)w Q (x), 
Q 

where {wq} is an orthonormal in L 2 (R) family of wavelets, such that wq is localized on 
Q. Define the Laplacian in the following way: 

a u = ]T2 2 ^)u q (;Vq(*), 

Q 

where 2~°^ is a side length of a dyadic cube Q. Katz and Pavlovic define the cascade 
operator as follows: 

(C(u,v)) Q = -2^UQVQ+2 Ej r- UQ V Q'> 

Q'ec 1 (Q) 

where Q is a dyadic cube with side length 2~i . The dyadic Navier-Stokes equation with 
hypo-dissipation is written as 

j t u + C(u,u) + v(A) a u = 0, 

where we include viscosity v, which is chosen to be one in 1141 . In terms of the wavelet 
coefficients uq, this equation can be written as 

jU Q {t)^-v2 2 ^u Q {t) + 2^u 2 Q {t)-2^ 11 u Q 

As it was proposed in [21 1 in the case v = 0, we simplify the model in the following 
way. Let Q\ be a dyadic cube with side length 2 _1 . Let v±(t) — UQ 1 (t) and vj(t) = 
UQj(t), where Qj is some dyadic cube in C^ 1 (Qi), j > 2. We will only consider the 
initial conditions for which uq(^o) = v j(to) f° r a U cubes Q e C^ 1 {Q\) for j > 2, and 
UQ(to) — for all dyadic cubes with side length larger than 2 _1 . Then, for every j > 2, 
we have UQ(t) = Vj(t) for all Q £ C 3_1 (Qi). Now, denoting vq — 0, we obtain the 
following system of equations for vj (t): 

j t v 3 {t) = -v2 2 ^v,(t) + 2 l irv]_ 1 - 2 3 2^v j v j+1 , j > 1, 
Finally, the change of variables 

u j (t)=2%v j (t/8). 
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reduces the equations to 

d 
It 

with v = f/8. 



— Uj = -v2 2aj Uj + 2 j uj_ 1 - 2 j+1 Uj u j+1 , j > 1, 



3. Functional setting 
Let us denote H = I 2 with the usual scalar product and norm: 

oo 

(U,V) := J~]unV n , \u\ := \J{u,u). 

71=1 

The norm \u\ will be called the energy norm. Let A : D(A) — > H be the Laplace operator 
defined by 

{Au) n = X 2an u n , n > 1, 
for some A > 1. The domain -D(-A) of this operator is a dense subset of H. Note that A is 
a positive definite operator whose eigenvalues are 

< X 2a < X 4a <X 6a <... 

Let H 1 = A^ 7 ^ 2a ^H endowed with the following scalar product and norm: 



{{u,v))j := ^2 X 2jn u n v n , \\u\\ 7 := ^/((u,u)) 7 . 

n=l 

In the special case 7 = a, let = H a = A~ X I 2 H and 

({u,v)) := ((u,u)) a , ||m|| := ||u|| a . 

This double norm ||w|| will be called the enstrophy norm. Note that we have an equivalent 
of the Poincare inequality 

H 2 <^\\u\\ 2 . 

Let also 

d s (u,v) := \u-v\, d w (u,v) := £ ^ J^-^L u ,v E H. 

n—l 1 1 

Here, d s is a strong distance, and d w is a weak distance that induces a weak topology on 
any bounded subset of H. Hence, a bounded sequence {u k } C H converges to u e H 
weakly, i.e., 

lim (u k , v) = (u, v), Vd e H, 

k^oo 

if and only if 

d w (u »0 as k — > 00. 

We also recall that if u fe — ► u weakly in as fc — > 00, then 

liminf \u k \ > \u\. 

Let 

C([0, T];H W ) := {«(•) : [0, T] -> u n (i) is continuous for all n} 
endowed with the distance 

dc([o,T];H v ){u,v) = sup d w (u(t) , v(t)) . 
te[o,T] 
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Let also 

C([0, oo); ff w ) :— {u(-) : [0, oo) — > H, u n {t) is continuous for all n} 
endowed with the distance 

^1 sup{d w (u(t),v(i)) :0<i<T} 

rfcdCoo);^) («, ») - ^ ^ 1+flup{dw(u(t)>7 , (t)):0 < t < r} • 

In this paper, the dyadic model of the Navier-Stokes equations will be written as 

(3 1) i ^u n + ^\ 2an u n - \ n ul_ 1 + \ n+1 u n u n+1 = g n , n = 1,2,3... 
[ uo = 0, 

for some parameter A > 1, the viscosity v > 0, the dissipation degree a > 0, and the force 
9 = (fflj 92) • • • )• F° r simplicity, assume that g is independent of time, g E H, and g n > 
for all n. 

For it = (u 1 ,U2, ■ ■ ■), the dyadic model can be written in a more condensed form as 

(3.2) —u + vAu + B(u, u) = g, 

at 

where 

-X n u n -iv n -i + X n+1 u n v n+ i, n = 2,3, . . . 
A 2 uiU2, n = 1. 

Clearly, the bilinear operator B enjoys the orthogonality property: 



(B(u,v)) n = 



OO 



(B(u,v),v) = 22 (-A n M„-i^n-i'y„ + X n+1 u n v n+1 v n ) 



= (-A ,l w„_iw„_iw„ + A"u„_iw n w„_i) 



= 0. 

Note that we always use a convention that uo = 0. 

Definition 3.1. A weak solution on [T, oo) (or (— oo, oo), if T = -oo) of ED is an 
H -valued function u(t) defined for t £ [T, oo), such that u n £ C 1 ([T, oo)) and u n (t) 
satisfies ( 13. 1 1 for all n. 

Note that since (B(u, u)) n has a finite number of terms, the notions of a weak solution 
and a classical solution (of a system of ODEs) coincide. Hence, the weak solutions will 
be often called solutions in the remainder of the paper. Note that if u(t) is a solution on 
[T, oo), then automatically u n 6 C°°([T, oo)). We say that a solution u(t) is strong (or 
regular) on some interval [T 1; T 2 ], if ||w(i)|| is bounded on [Ti, T?\. A solution is strong on 
[Ti, oo), if it is strong on every interval [Ti,T 2 ], T 2 >T X . 

Definition 3.2. A Leray-Hopf solution of ( 13. U on the interval [T, oo) is a weak solution 
of d3 . 1 1 on [T, oo) satisfying the energy inequality 

(3.3) \u(t)\ 2 + 2is[ \\u{T)\\ 2 dT <\u(t )\ 2 + 2 [ (g,u(T))dT 

for all T < to < t, to a.e. in [T, oo). The set Ex on which the energy inequality does not 
hold will be called the exceptional set. 
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Note that the complement of the exceptional set Ex coincides with the set of points 
of strong continuity from the right. Later we will prove that every solution u(i) with 
u n (T) > is a Leray-Hopf solution on [T, oo), and that the energy inequality for such a 
solution is satisfied starting from any time to > T, i.e., Ex = 0. 



4. A PRIORI ESTIMATES AND EXISTENCE OF WEAK AND STRONG SOLUTIONS 

We start with some a priori estimates. 
Energy estimates. Formally taking a scalar product of the equation J3 . 1 1 with u, we obtain 

lj t \u\ 2 <-v\\u\\ 2 + \g\\u\ 

<- V \ u \ +_ H +_ 

= -"-\u\ 2 + K 
V 1 2v 



Using Gronwall's inequality, we conclude that 



(4.1) | u ( t )|2< e -*| u (o)| 2 + %(l-e 



2 



-vt\ 



V 



Hence, B = {u E H : \u\< R} is an absorbing ball for the Leray-Hopf solutions, where 
R is any number larger that \g\jv. Note that this result will later follow rigorously from 
the energy inequality. 

Next, taking a limit of the Galerkin approximation, we will prove the existence of Leray- 
Hopf solutions to i3.l\ . 



Theorem 4.1. For every u° £ H and g € H, there exists a solution of (13 . 11 with u(0) = 
u . Moreover, the energy inequality 

\u{t)\ 2 + 2v [ \\u(r)\\ 2 dT<\u(t )\ 2 +2 f\g,u(T))dT 

J to Jto 

holds for all < to < t, to a.e. in [0, oo). 

Proof. Let u° 6 H. We will show the existence of a weak solution by taking a limit of 
the Galerkin approximation u k (t) = (u\(t), . . . , u\{t), 0,0,...) with it„(0) = u° for 
n = 1, 2, k, which satisfies 



(4.2) 



-«* + v\ 2an u\ - X n (u k n _ 1 ) 2 + X n+1 u k n u k n+1 = g n , n <k — l, 



where it§ = 0. First, note that the energy estimate ( 14. 1> obviously holds for u k (t). Hence, 
from the theory of ordinary differential equations we know that there exists a unique solu- 
tion u k (t) to (14. 2> on [0, oo). Next, we will show that a sequence of the Galerkin approxi- 
mations {u k } is weakly equicontinuous. Indeed, thanks to the energy estimate ( 14. 1> . there 
exists M, such that 

Un{t) < M, Vn, k, t > 0. 
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Therefore, 

\u k n (t) - u k n (s)\ < J* (-v\ 2an u k + A>*_i) 2 - + 9n) dr 

< {v\ 2an M + \ n M 2 + X n+1 M 2 + g n )\t - s\, 
for all n, k, t > 0, s > 0. Thus, 

<®-<(s)\ 



„_,- - + !<(*)-<(*) I 

<c\t-s\, 

for some constant c independent of k. Hence, {u k } is an equicontinuous sequence of func- 
tions in C([0, oo); H w ) with bounded initial data. Therefore, the Ascoli-Arzela theorem 
implies that {u k } is relatively compact in C([0, T]; H w ) for every T > 0. By a diagonal- 
ization process it follows that {u k } is relatively compact in C([0, oo); H w ). Hence, passing 
to a subsequence, we obtain that there exists a weakly continuous H -valued function u(t), 
such that 

(4.3) u kj ^u as kj — > oo in C([0, oo); H w ). 

In particular, uf/ (t) — > u n (i) as fcj — > oo, for all n, t > 0. Thus, it(0) = u°. In addition, 
note that 

i# (*) = t# (0) + { {-vX^vft + \ n (v%_ 1 ) 2 - A» +1 l#t# +1 + g n ) dr, 
Jo 

for n < kj — 1. Taking the limit as kj — ► oo, we obtain 

u»(*)=«n(0) + / (-^A 2m V„ + \ n ul_ t - X n+1 u n u n+1 + g n ) dr. 



Since u n (t) is continuous, it follows that u n <= C 1 ([0, oo)) and satisfies J3. II . 

It remains to prove that u(i) satisfies the energy inequality. Note that u kj (t) satisfies 
the energy equality 



\u^{t)\ z + 2v ||u^(T)|rdr=|tt^(to)r + 2 / (g,u^(r))dT, 
Jto J t 

for all t > to > 0. Hence, the sequence {u kj } is bounded in L 2 ([to, t]; V) for all t > to > 
0. This, together with (14.31 . implies that 



\u ki {r) ~ u(t)\ 2 dr — > 0, as fej — > oo, 

for all i > <o > 0. In particular, |it fej (i)| — > |it(i)| as fcj — > oo a.e. in [0, oo). Take any 
< > for which |u ■'(io)l ~* l u (^o)l as &j ~~ * °°- For every N > 0, we have 

|^(t)| 2 +2z, /' £ \ 2an u k ;{r) 2 dr < |^(t )| 2 +2 f (g,u k > (t)) dr. 
Jto n < N Jto 

Since u fej (t) — > u(t) weakly in H as fej — * oo for all time t > 0, we have that 

| U (t)| 2 + 2^ f X 2an u n (r) 2 dr<\u(t )\ 2 + 2 f (g,u(r)) dr. 
Jt „</v Jto 
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Finally, taking the limit as N — > oo and using Levi's theorem, we obtain 



\u(t)\ 2 + 2v / \\u{T)\\ 2 dr< Ki )| 2 + 2 / (g,u(r))dT, 



it, 



In 



for all < to < t, to a.e. in [0, oo). 



□ 



Note that this was a classical proof from the theory of the NSE. Using the fact that there 
is no backward energy transfer, we can actually show that every solution with u n (0) > 
is a Leray-Hopf solution and, moreover, is continuous from the right in H for all time. 

Theorem 4.2. Let u(t) be a solution of (13.11 with it„(0) > 0. Then u n (t) > for all 
t > 0, and u{t) satisfies the energy inequality 



(4.4) 



\u{t)\ 2 + 2v / \\u(T)\\ 2 dr< Kt )| 2 + 2 / (g,u(r))dT, 



In 



In 



for all < t < t. 

Proof. A general solution for u n (t) can be written as 

(4.5) Un(t) = u n (0) exp (- f v\ 2an + X n+1 u n+1 (r) dr 



+ / exp 



uX Zan + A "+l Wn+l(r ) dr ^ + AX-l(*)) (fe. 



Recall that g n > for all n. Since u n (0) > for all n, then u n (t) > for all n, t > 0. 
Hence, multiplying (13 . 1 1 by u n , taking a sum from 1 to N, and integrating between to and 
t, we obtain 

N N „t N 

5>„(t) 2 - ]>>„(to) 2 + 2z, / ^A 2ttn Mll (r) 2 dr 

i 1 «/ tn 1 



71=1 



= -2 / A Ar+1 u 2 v itiv+i + 2 / g n u n dr 

Jto J to n= i 

,t N 

- 2 / X! gnUn dr - 

Jtn i 



Taking the limit as N — > oo, we obtain (I4.4i . 

Enstrophy estimates. We obtain the following estimate for the nonlinear term: 



□ 



\(B(u,u),Au)\ 



oo 

j~^2a(n+l 



)-A 



2an 



A u n u n+ i 



= (A" - A -a ) 



^ A^ +1 ^ 2 A Q(n+1) u„+i 



< 



c b (max|A a "u„|) V X^ a+1 ^ n u 2 n 

n. i J 



< c h \\u\\J2 X^ +1 



)n u 2 
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where Cb = A" — A a > 0. When a <E [1/3, 1], Holder's inequality implies 

l/a-l| 4 l/2 i3-l/a 



\(B(u,u),Au)\ < c^ulWAul 1 ^- 1 ^ 



= c^Au^l^WuW 4 - 1 / 01 . 

Choosing u to have only two consecutive nonzero terms, it is easy to check that these 
estimates are sharp. Moreover, when a — 2/5, we have 

\(B(u,u),Au)\ < c b |A u | 3 / 2 |M| 3 / 2 , 

which is the same as the Sobolev estimate for the inertial term of the 3D NSE (see, e.g., 
)6|l20l '). Therefore, taking a scalar product of the equation J3 . 1 1 with Au and using Young's 
inequality, we obtain 

\j t \\uf < -v\Au\ 2 + c h \Au\^ 2 \\uf' 2 + (g, Au) 

< - v \Au? + %\Au\* + |^N| 6 + A|. 9 | 2 + ^ A u\ 2 

a Riccati-type inequality for ||u|| 2 . Hence, the model has the same enstrophy estimate as 
the 3D NSE, similar properties, and the same open question concerning the regularity of 
the solutions in the case a = 2/5. 

Another interesting case is a = 1/3. Then we have 

(4.6) \{B(u,u),Au)\<c h \Au\ 2 \\u\\, 

which corresponds to the 4D Navier-Stokes equations. 

Now consider the case where a > 1/3. Formally, we have 

\j t \M 2 < -v\M 2 + Cbl^l^lMI 4 - 1 /" + (g,Au) 

< ^\ Au \ 2 + ^\Au\ 2 + c\\u\\^ + ^\g\ 2 + "-\Au\ 2 

V o 8a-2 3 o 

= --| J 4 u | 2 + c || u ||^T + -| ff | 2 , 

for some constant c > 0. This means that if the initial data is in V, then u(t) remains 
bounded in V for some time T. Applying the above estimate to the Galerkin approximation 
and taking a limit, we immediately obtain the following local regularity result. 

Theorem 4.3. If a > 1/3, then for any u° G V there exists a strong solution u(t) to ( 13. H 
on some time interval [0, T], T > with u(0) = u a . 

Finally, consider the case a > 1/2. In this case the enstrophy estimate implies 

(B(u,u),Au) < c h \Au\\\u\\ 2 . 

Therefore, formally, we have 

lj t \\u\\ 2 < -v\Au\ 2 + c h \Au\\\u\\ 2 + (g,Au) 

< _ v{Au \2 + ^ Au \ 2 + M|M| 4 + A|.g| 2 + V -\Au\ 2 
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This is again a Ricatti-type inequality. Assume that u(t) is a strong solution on some 
interval (0,**), and ||u(£)|| -> +00 as t -> t*-. Then 

||u(*)|| 2 > j^r t , 0<t<t*, 

for some positive constant c. However, this means that ||it(£)|| 2 is not locally integrable, 
which is in contradiction with the energy inequality. Hence, if the initial data u° e V, then 
||u(f)|| is bounded on every interval [0, T], T > 0, and we have the following. 

Theorem 4.4. If a > 1/2, then for any u° (£ V there exists a strong solution u(t) to ( 13. H 
on [0, 00) with u(0) = u . 

5. Blow-up in hnite time 

Let a < 1/3 and 7 > 0. In this section we will prove that every solution u(t) with 
large enough ||u(0)|| 7 blows up in finite time in H 1 ^^ 1 norm. The idea is the following. 
Taking a scalar product of the equation with A^^u, we obtain 

~IM£ = -HMI« +7 + «), ^ 7/Q «) + (s, 

In order to show a blow-up, we, in some sense, will invert the Sobolev estimates for the 
nonlinear term. Note that 

(B(u,u),A^ a u) ~Y. Xil+2l)nu2 n u ^- 

n 

If u n > is monotonically decreasing in n, then 

Obviously, this is not true in general. For example, if u n — for even n, then (B(u,u),Ai/ a u) 
0. However, we will prove that a similar estimate holds if we use the following function 
instead of _ff 7 -norm: 

H(t) :=\\u(t)\\ 2 1 + cJ2^~< n (u n u n+1 )(t), 

n 

for some constanat c > 0. More precisely, we will show that if ||u|| 7 is large enough, then 
1 d 



, di H>-v\\u\\l + , + j:^ )n ul 



provided that a < 1/3. We will start with the following estimate. 

Lemma 5.1. If a < 1/3, then for any 7 G (0, 1 — 3a) there exists a positive constant A, 
such that 

71=1 
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Proof. Let e := 2 - 6a - 27 > 0. Note that A~ e < 1. Let 
Holder's inequality with p = 3 and g = 3/2 implies 



71=1 



< 



1/3 



2/3 



£ A (1+2 7 )„|, 
\n—l 

\ 2 /3 



Hence, 



= A ~ 2/3 E A(1_ 

\n=l / 
00 

£A^> n | 3 >AH 3 a+71 



n=l 



which concludes the proof. 



□ 



Lemma 5.2. Let u(t) be a solution to ( 13. \\ and 7 > 0. Let be continuous on 

[0, 00). Then the function 



(5.1) 



n=l 



is continuous on [0, 00). 

Proof. Let v n (t) — X 2ryn (u n u n +i)(t). First, note that due to the Cauchy-Schwarz in- 
equality, the function i5.li is less than or equal to and, consequently, is bounded 
on every interval [a, b], < a < b. Let to > 0. Since is continuous at t = to, it 
follows that 

lim limsup V \ 2ln u 2 n {t) = 0. 

l ^ t0 n=N 



Therefore, 



lim limsup > v n (t) — 0, 



-N 



which means that ( 15. It is continuous at t = to. Indeed, since v n (t) is continuous for every 
71, we have 



lim sup 

t->t 



£««(<) - £«n(*o) 



n=l 



lim lim sup 



N~l 



N-l 



-N 



n—1 n—1 n—N 

Similarly, the continuity of (15. 1> from the right holds at i = 0. 
Now we proceed to our main result. 



□ 
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Theorem 5.3. Let u(t) be a solution to ( 13. H with u n (0) > and a < 1/3. Then for every 
7 > 0, there exists a constant M("f), such that \\u(t)\\^ 3+ ^ is not locally integrable on 
[0, oo), provided ||u(0)|| 7 > M(-j). 

Proof. Since ||u|| 7l < |M| 72 for 71 < 72, it is enough to prove the theorem in the case 
< 7 < min{l/3, 1 — 3a}. Given such 7, let u(t) be a solution to (13. 1> . such that 
ll w (*)ll 1/3+7 i s integrable on [0, T] foreveryT > 0. We will show that ||it(0)|| 7 is bounded 
from above by a constant dependent of 7. 

Note that u n {t) > for all n, t > due to Theorem l4.2l First, we obtain 

J n =l J ° n=l 

Jo \n=l 



(5.2) 



\ 3/2 

<2 I IVAl^V dr 



<2 / \\u(t)\\ 3 1/3+ dT 







< 00, 



for all T > 0. Thus, 



£ A( 1+2 ^(u^ n+1 ) (t) and £ X^) n u n {tf 

n—1 n—1 

are locally integrable on [0, oo). In addition, since a < 1/3, we have 

|N*)lla+ 7 < IN*)lll/3+r 

which implies that ||w(i)||a +7 is locally integrable on [0, oo). 

Now note that if u n < ^u n+ i, then u n u 2 l+1 < \u^ l+l . Otherwise, u n u^ +1 < 
2u 2 l u n +i- Hence, 

(5.3) u n u 2 n+l < \u 3 n+1 + 2u I 2 l M n+ i, n € N. 

This also implies that 

U n U n+ iU n+2 < \ul t U n+1 + \u n+1 U^ +2 

* ^ 12 13 2 

< 2 U n u n+l + J u n+2 + u n+l u n+2i 

for all neN. 

From ( 13. 1> we obtain 

^K^+i) = - ^(A 2mi + X 2a ^)u n u n+1 
at 

+ X n u 2 n ^ 1 u n+1 - \ n+1 u n u 2 n+1 

+ \ n+1 ul - \ n+2 U n U n+1 U n+2 

+ g n u n+ i + g n +iu„. 
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This, together with inequalities J5.3I and J5.4I . implies that 

j t {u n u n+1 ) + v(l + X 2a )X 2an u n u n+1 

+ 2X n+1 u 2 n u n+1 + \\ n+2 u 2 n u n+1 + \ n+2 u 2 n+l u n+2 

> \n+l 3 lxn+1 3 1 xn+2 3 

— " U n 2 A U n+1 ~ 1 A U n+2- 

Multiplying it by A 27 ™, taking a sum from 1 to oo, and integrating between and f, we get 

oc oo 

J2 A 27 "K Mn+1 )(t) - ]T A 27 "K«„+i)(0) 

n— 1 Ti—l 

pt 

+ is{l + \ 2a ) / ^ (a+l)n UnU n+1 dT 

r t oo 

(5.5) +(2A+iA 2 + A 1 - 27 ) / VA^+^V^n+irfr 

> (A-iA- 27 -iA- 47 ) / ^A( 1+2 ^ dr 

•/O n=l 

•'° n=l 

for all £ > 0. On the other hand, we have the following equality for the nonlinear term: 

oo oo 

-(B(u, u),A^ a u) = \ 2ln+n u 2 n _ lUn - \ 2ln+n+1 u 2 n u n+1 

71—1 71 — 1 

OO 



n=l 

where c\ — A 27+1 — A > 0. Now, multiplying 13. II by A 27 ™w„, taking a sum from 1 to 
oo, and integrating between and t, we obtain 

(5.6) \\u(t)\\ 2 -\\u(0)\\ 2 +2is / || U (r)|| 2 +7 dr 

Jo 

t oo 



pt ^o pt ^o 

2 Cl / V \^ +2 ^ n u 2 n u n+1 dr + 2 V A 27 ™ 5 „ U „ dr, 
Jo __i Jo i 



for i > 0. Note that the term with the force is integrable because j < 1/3. In particular, 
(15.61 yields that || it (i)|| 2 is continuous on [0, oo). Denote 

oo 

H(t) := || U (t)|| 2 + c 2 \ 2 ~< n (u n u n+1 )(t), 

n=l 

where c 2 = 2ci/(2A + A 2 /2 + A 1-27 ). Thanks to Lemma |5~2*1 H(t) is continuous on 
[0, oo). We will show that H(t) is a Lyapunov function, i.e., H(t) is always increasing. 
Moreover, we will see that H(t) blows up in finite time. Indeed, multiplying ( 15.51 by c 2 
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and adding (15.61 . we get 

pt ft °° 

H(t)-H(Q)>-2u \\u{r)\\l +1 dT-uc s V X 2( - a+ ^ n u n u n+1 dr 
Jo Jo n=1 

J ° n=l 

where C3 = (1 + A 2 ")c2. Due to Lemma lBTl there exists a constant A > 0, such that 



n=l 

In addition, the Cauchy-Schwarz inequality implies 



\ 2{a+l)n u n u n+1 < \\u\ 



2 

a+7' 



Therefore, we obtain 

(5.7) H(t) - H(0) > -v{2 + c 3 ) £ ||u(r)||2 +7 dr + ^ J* \Hr)f a+1 dr, 
for t > 0. Note that 

||«(t)||* <ff(t)<(l + 02)||u(t)||*. 

In particular, 



|K*)||c« +7 > 

Let 



H{t) 
l + c 2 ' 



MM := MHC3) ^ . 
AAc 2 

Assume that H(t) > M 2 on [0,t] for some £ > 0. Then we have that ||w(t)|| q+7 > 
8u(2 + c^)/{A\c2) on [0, t] and, consequently, ( 15. 7> yields 

(5.8) - #(0) >c f H{rf/ 2 dr, 

Jo 

where c = AAc 2 /8. Now assume that ||u(0)|| 7 > M. Then H(0) > M 2 and Ol holds 
for some small time t > 0. Then (15. 8> automatically holds for every t > 0. 

Note that J5.8i is a Riccati-type inequality. It is easy to see that H(t) blows up in finite 
time. Indeed, let y(t) be the solution to the Riccati equation 

y'(t) = cy(tf/ 2 , y(0) = |H(0). 

Then for some t* > 0, we have that y(t) — > 00 as £ — > £* — . Consider 

w(t)=H(t)-y(t). 

It is easy to check that the function satisfies the following integral inequality: 

w(t) - w(0) >cf w(t) 3/2 dr, 
Jo 

for all t > 0, such that w(t) > on [0,t]. Note that w(0) > and u>(t) is continuous. 
Thus, w(t) > for all t € [0, t*). 

Now, since y(i) blows up in finite time, H(t) also blows up in finite time, which con- 
tradicts the fact that H(t) is continuous on [0, 00). Hence, ||u(0)|| 7 < M. □ 
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6. Global attractor 

In Section|4]we showed that the dyadic model possesses an absorbing ball with a radius 
R larger than \g\/u. Let X be a closed absorbing ball. 

X := {u e H : \u\ < R}, 

which is compact in d w -metric. Then for any bounded set K C H, there exists a time to, 
such that 

u(t) el, Vf > t , 

for every Leray-Hopf solution u(t) to ( 13. 11 with the initial data u(0) G K. 

We recall the definition of an evolutionary system £ from [4| (see also |3|). Let 

T :={I : I = [T, oo) C R, or / = (-oo, oo)}, 

and for each I C T let denote the set of all X-valued functions on I. A map £ that 
associates to each / 6 T a subset £(I) C T will be called an evolutionary system if the 
following conditions are satisfied: 

(1) £([O,oc))^0. 

(2) £(I + s) = {u(-) : u(- - s) E £(I)} for all s E M. 

(3) {u(-)|/ 2 : u{-) € £(h)} C £{h) for all pairs of h,I 2 € a such that J 2 C h. 

(4) £({-oo, oo)) = {u(-) : u(-)|[t,oo) € £([T, oo)) VT e M}. 

Let 

:= {u(t) : u(0) eA,u£ £{[0, oo))}, 

R(t)A := {u(t) : u(0) £A,ue £((-oo, oo))}, A c X, t > 0. 

For A C X and r > 0, denote r) = {u : d m (A, u) < r}, where • = s, w. Now 

we define an attracting set and a global attractor as follows. 

Definition 6.1. A set A C X is a d. -attracting set (• = s, w) if it uniformly attracts X in 
d.-metric, i.e., for any e > there exists t , such that 

R{t)X c B,(A,e), Vt >t . 

A set A C X is invariant if R(t)A = A for all t > 0. A set A, C X is a d.-global 
attractor if A, is a minimal d. -closed d. -attracting set. 

The following result was proved in |4|: 

Theorem 6.2. The evolutionary system £ always possesses a weak global attractor A w . 
In addition, if £([0, oo)) is compact in C([0, oo); H w ), then 

(a) A w = {u° : u = u(0)for some u € £((— oo, oo))}. 

(b) A w is the maximal invariant set. 

For the dyadic model, we define £ in the following way. 

£([T,oo)) :—{u : u(-) is a Leray-Hopf solution on [T, oo) 
and u(t) e XVt e [T,oo)}, Tel, 
£ ((oo, oo)) :={u : u(-) is a Leray-Hopf solution on (— oo, oo) 
and u(t) elVie (-oo, oo)}, 
where X is the phase space defined in the beginning of the section. Clearly, £ satisfies 
properties (l)-(4). Then Theorem 16 . 21 immediately yields that the weak global attractor 
„4 W exists. In order to infer that A w is the maximal invariant set, we need the following 
result. 
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Lemma 6.3. £ ([0, oo)) is compact in C([0, oo); H w ). 
Proof. Take any sequence u G £([0, oo)). First, note that 

u k n (t)<R, Vn,k,t>0. 

Therefore, 

\u k n {t) - u k n {s)\ < (v\ 2an R + \ n R 2 + \ n+1 R 2 +g n )\t - s\, 
for all n, k, t > 0, s > 0. Thus, 

^2(« 2 ) i + |<(t)-<( s )| - 

for some constant c independent of k. Hence, {u k } is an equicontinuous sequence of 
functions in C([0, oo); H w ) with bounded initial data. Therefore, Ascoli-Arzela theorem 
implies that {u k } is relatively compact in C([0, T]; -ff w ) for all T > 0. Using a diagonal- 
ization process, we obtain that {u k } is relatively compact in C([0, oo);H w ). Hence, there 
exists a weakly continuous i/-valued function w(i) on [0, oo), such that 

(6.1) u kj — > u as fej- — ► oo in C([0, oo); H" w ), 

for some subsequence fcj . In particular, 

|it(t)| < liminf \u k i{t)\ < R, t>0, 

kj — »00 

i.e., u(t) G X for all t > 0. 

In addition, since M fej (t) is a solution to ( 13. U . we have 

t# (t) = Uj? (0) + \\-v\ 2an U k - + X n (u k n U) 2 - ^ n+lu n u n+l + 9n) *T, 

Jo 

for all n. Taking the limit as fcj — * oo, we obtain 

u n {t) = u n (0) + f (-is\ 2an u n + X n u 2 n _ 1 - X n+1 u n u n+1 + g n ) dr, 
Jo 

for all n. Since u n (t) is continuous, u n G C 1 ([0, oo)) and satisfies J3 .11 . 

In order to infer that u G £{[0, oo)), it remains to prove that u(t) satisfies the energy 
inequality. Note that \u kj — > |«(t)| as fcj — > oo a.e. in [0, oo). Since u kj G £([0, oo)), 
it satisfies the energy inequality starting from any to that is not in the exceptional set of 
measure zero. Let Ex be the union of the exceptional sets for all u kj . Note that Ex is of 
measure zero. Take any to ^ Ex for which \vr* (t ) \ — > |u(to)| as &j — > oo. Then 

|^(t)| 2 + 2^ /' t ||^(T)|| 2 dr<|^(to)| 2 + 2 f (g,u k > (t)) dr, 

J to J to 

for all t > t n . Hence, 

|w fe ^)| 2 + 2i/ / ^ A 2Q "u^ (t) 2 dr < \u k i{t Q )\ 2 + 2 f\g,u k >(T))dT. 

J to „< Ar J to 

Since u fej (t) — > u(t) weakly in as fej — * oo for all time i > 0, we have that 

\u{t)\ 2 + 2v f Y, X 2an u n (r) 2 dT<\u(t Q )\ 2 + 2 f\g,u(r))dT. 

J t n<N J t 
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Finally, taking the limit as N — > oo and using Levi's convergence theorem, we obtain 



\u(t)\' + 2v / \\u(T)\\ 2 dr< \u(t )\ 2 + 2 / (g,u(T))dT, 

J t J to 

for all < to < t, to a.e. in [0, oo). Hence, u € £([0, oo)), which concludes the proof. □ 

Now Theorem 16.21 implies that the weak global attractor A w is the maximal invari- 
ant set that consists of the points that belong to complete trajectories, i.e., trajectories in 
£({oo, oo)). Moreover, using ( 14. 5\ . one can show that u n > for every u S A w . Consider 
now the case a < 1/3. It is easy to show that for every 7 £ (0, 1 — 3a), we can take g\ 
large enough, so that for every solution u(t) and every t > 0, we have \u(t)\ > M{^) for 
some t € [t, t+ 1]. Thanks to Theorem l5.3l this means that A w is not bounded in i/ 1 / 3+7 . 
It is an open question whether A w is bounded in V. 

We will now proceed to study the question whether A w is also a strong global attractor. 

Theorem 6.4. Let a > 1/2. 77ien every Leray-Hopf solution u(t) of ( 13. H satisfies the 
energy equality 



\u{t)\ 2 + 2 V \\u(r)\\ z dT = \u(t a )\' + 2 {g,u(T))dT, < t a < t 



to 



Proof. Let u(t) be a Leray-Hopf solution of J3. II . Thanks to the energy inequality ( 14. 41 . 
\\u(t)\\ 2 is locally integrable. Then we obtain 

^ t o n=l n=l 

ft 00 



<2 sup \u(s)\ / ^A"<dr 
se[to.t] J to n =l 

< 2R f \\u(t)\\ 2 dr 

Jtn 



< OO, 



for < to < t. Hence, 

/ \ H+1 u 2 n u n+ i dr — > 0, as n — > 00. 

J to 

Multiplying (13- II by u n , taking a sum from 1 to N, and integrating between to an d t, we 
obtain 



N N .t N 

n=l 71=1 •' t ° 71=1 



-2 / A^u^ttjv+i dT + 2 
•/to 



7t " n=l 



Finally, taking the limit as AT — * 00, we arrive at 

\u(t)\ 2 + 2v[ \\u(t)\\ 2 dr = \u(t Q )\ 2 + 2 I ( 5)U (r))dr, < i < t. 



to 



□ 
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In 1 4 1 it was proved that the asymptotic compactness of the dynamical system £ implies 
that the strong global attractor A s exists, is strongly compact, and coincides with A w - In 
the case where the evolutionary system consists of the Leray-Hopf weak solutions to the 
3D NSE, the continuity of the complete trajectories, i.e. trajectories on A w , implies the 
asymptotic compactness of £ (see also 1 1 1 and 1 19 1 for similar results). In 1 3 1 this result was 
proved for an abstract evolutionary system satisfying the energy inequality. It immediately 
implies the following. 

Corollary 6.5. Let a > 1/2. Then A w is a strongly compact strong global attractor. 

Note that if a > 1/2, then, thanks to Theorem l4.4l for every initial datum in H there 
exists a regular solution on [0, oo). Moreover, it can be shown that such a solution is unique 
in the class of all Leray-Hopf solutions. Hence, Corollarv l6.5l can also be obtained using a 
classical theory of semiflows. It is an open question whether the continuity of the complete 
trajectories and, consequently, the existence of the strong compact global attractor holds 
for a < 1/2. 

ACKNOWLEDGMENT 

The author would like to thank Susan Friedlander and anonymous referee for helpful 
comments and suggestions. 

References 

[1] J. M. Ball, Continuity properties and global attractors of generalized semiflows and the Navier-Stokes equa- 
tions, J. Nonlinear Sci. 7 (1997), 475502. Erratum: J. Nonlinear Sci. 8 (1998), 233. 
[2] L. Biferale, Shell models of energy cascade in turbulence, Annu. Rev. Fluid Mech. 35 (2003), 441468. 
[3] A. Cheskidov, Global attractors of evolutionary systems, preprint, (2006). 

[4] A. Cheskidov and C. Foias, On global attractors of the 3D Navier-Stokes equations, preprint, (2006). 

[5] A. Cheskidov, S. Friedlander, and N. Pavlovic, The global attractor of an inviscid dyadic model for the fluid 

equations, preprint, (2006). 
[6] P. Constantin and C. Foias, Navier-Stokes Equation, The University of Chicago Press, Chicago, 1989. 
[7] P. Constantin, B. Levant, and E. Titi, Analytic study of the shell model of turbulence, preprint, (2005). 
[8] E. I. Dinaburg and Ya. G. Sinai, A quasi-linear approximation of three-dimensional NavierStokes system, 

Moscow Math. J. 1 (2001), 381-388. 
[9] S. Friedlander, and N. Pavlovic, Blow up in a three-dimensional vector model for the Euler equations, 
Comm. PureAppl. Math. 57 (2004), 705-725. 
[10] S. Friedlander, and N. Pavlovic, Remarks concerning modified Navier-Stokes equations, Discrete and Con- 
tinuous Dynamical Systems 10 (2004), 269-288. 
[11] U. Frisch, Turbulence. The Legacy of A. N. Kolmogorov, Cambridge University Press, Cambridge, 1995. 
[12] E. B. Gledzer, System of hydrodynamic type admitting two quadratic integrals of motion, Soviet Phys. Dokl. 
18 (1973), 216-217. 

[13] N. H. Katz and N. Pavlovic, A cheap Caffarelli-Kohn-Nirenberg inequality for the NavierStokes equation 

with hyper-dissipation, GAFA 12 (2002), 355-379. 
[14] N. H. Katz and N. Pavlovic, Finite time blow-up for a dyadic model of the Euler equations, AMS Tran. S 

0002-9947(04)03532-9 (2004). 
[15] A. Kiselev and A. Zlatos, On discrete models of the Euler equation, IMRN 38 (2005), 2315-2339. 
[16] V. S. Lvov, E. Podivilov, A. Pomyalov, I. Procaccia, and D. Vandembroucq, Improved shell model of 

turbulence, Phys. Rev. E (3) 58 (1998), 1811-1822. 
[17] A. M. Obukhov, Some general properties of equations describing the dynamics of the atmosphere, hy. Akad. 

Nauk SSSR Ser. Fiz. Atmosfer. i Okeana 7 (1971), 695-704. 
[18] K. Ohkitani and M. Yamada, Temporal intermittency in the energy cascade process and local Lyapunov 

analysis in fully-developed model of turbulence, Progr. Theoret. Phys. 81 (1989), 329-341. 
[19] R. M. S. Rosa, Asymptotic regularity condition for the strong convergence towards weak limit sets and weak 

attractors of the 3D Navier-Stokes equations, preprint. 
[20] R. Temam, Navier-Stokes Equations: Theory and numerical analysis, North-Holland, Amsterdam, 1984. 
[21] F. Waleffe, On some dyadic models of the Euler equations, AMS Proc. to appear. 



BLOW-UP FOR THE DYADIC MODEL 



University of Michigan, Department of Mathematics, Ann Arbor, MI 48 109 
E-mail address: acheskid@umich . edu 



